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Abstract. Monodisperse circular disks have been collectively packed in confined 
geometries using a Monte Carlo method where the compaction is propelled by two- 
dimensional stochastic agitation. We have found that confinement (i.e., finite-size plus 
surface effects) determines the symmetry of the packed configurations together with 
the size evolution of the probability density of the packing fraction. For the particular 
case of small systems in square containers, the probability density of the packing 
fraction shows several well-defined peaks, depending on the system size, for both hard- 
wall and periodic boundary conditions. High-symmetry configurations (other than the 
nxn square arrays) with non- negligible occurrence probabilities are found as a direct 
consequence of monodispersity and confinement. 
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1. Introduction 

Confined monodisperse systems have recently attracted significant attention for their 
potential in technological applications and their role in explaining fundamental questions 
on self organization p~H9]. For instance, current experimental techniques are able to 
produce extremely monodisperse foams, bubbles, and colloidal nanoparticle suspensions 
that are stable against coalescence and coarsening in a time span of several days. This 
renaissance in monodisperse systems has been spurred by interest in manipulating fluids 
in micrometric channels (microfluidics) pp. The possible applications of microfluidics 
range from chemical and biological analysis to optics to information technologies [2]. 
In addition, there are striking analogies between monodisperse bubbles and granular 
materials that extend to both jamming and flow [3]. 



Monodisperse systems tend to self organize in nontrivial ways under confinement pE]. 
Notably, monodisperse colloidal particles have shown a rich configurational landscape in 
two-dimensional channels [5116] and planar patterning on surfaces [7HH] . The possibility of 
controlling the state of order in monodisperse particles by confinement has applications 
ranging from photonic crystals to drug-delivery to sensor technologies [4] . 

In the cases discussed above, different mechanisms were employed to activate and 
sustain the self-assembly process. For example, a laminar flow was used to stabilize 
the novel crystal structures of colloidal particles in microchannels in Ref. [6], whereas 
an electrostatic field aided the self-assembly of electrodeposited particles of Ref. [9]. A 
very attractive and viable process for self-assembly is vibration, either in the form of 
thermal or mechanical agitation. The tapping experiments in granular materials are 
well-known examples of the latter [10]. Vibration is at the core of the self-assembly 
of micro-to-millimeter objects with a variety of functions from computing devices [H] 
to semiconducting devices (LEDs) [12] to three-dimensional electric circuits [13]. Quite 
recently, a high- yield vibration-based method to assemble magnetically interacting cubes 
has been reported [14J. 

Here we propose, and analyze the consequences of, a simple model for the self- 
assembly of two-dimensional confined hard disks under stochastic vibration (Sec. |2]). In 
contrast to previous works on confined systems of monodisperse disks or spheres, which 
have mostly been concerned with the random close packing limit (see, e.g., Ref. [15] for a 
recent and excellent survey of the relevant literature), our goal here is to investigate the 
impact of the boundary conditions on the distribution of configuration states. Special 
emphasis is given to identifying nontrivial high-symmetry configurations and their yield 
rates — a subject of interest for applications on devices whose functionality is based upon 
different states of order jUEHH] (Sec. |3]). 

2. Model approach and compaction protocol 



We consider the following scenario: Two-dimensional monodisperse disks are distributed 
randomly, following a uniform distribution, in a square container with hard walls. 
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The density of the system (initially diluted) is uniformly increased by compressing 
the container (or equivalently by expanding the disks) with a slow yet monotonic 
rate 7. During this compaction process, the particles are continuously agitated in the 
horizontal plane with frequency u and amplitude Aj. Each particle % is individually and 
randomly shaken. The shaking amplitude is random, following a Gaussian distribution 
characterized by zero mean and standard deviation 07. The direction of the agitation 
is also random but, in contrast to the amplitude, the shaking direction follows a flat 
random angular distribution over 2ir. These choices are in place to guarantee that the 
particles sample different positions within the container with no preferential direction, 
i.e., a truly two-dimensional stochastic sampling. The compaction process continues, 
resolving particles overlaps through the stochastic agitation, until no further increase 
of the density is possible. A protocol based on this approach consistently appears to 
unbiasedly survey the configurational space, retrieving rare events with small, yet finite 
probability (see below). 

We have implemented our model in a Monte Carlo approach where the particles 
are frictionless and represented by regular polygon^. We assume a frictionless hard- 
body interaction between the particles and the walls of the container. The frequency of 
vibration u becomes the number of shakes per particle. We have scanned the parameter 
space for the least-biased set of values for the growth rate 7, shaking frequency u, and 
&i [IE] — note that least-biased in this context means the subset of values where the error 
bars were minimized for at least 5 runs, and the packing fraction had reached a steady- 
state value with respect to the variable. We have found that a compaction protocol with 
7=10~ 5 , Co>=1000, and ai/d=l (d is the disk diameter), offers a reasonable sampling of 
the configurational space without favoring particular disks arrangements. 

The significant length scales are the diameter of the disks and the size of the 
container. In an inflationary packing scheme, as the one used here, the ratio of those 
length scales becomes equivalent to the total number of particles (see, for example, 
Ref. |17j). In the following, we investigate how n 2 (=9-1024) disks are packed in square 
containers. Doing so allows us to pinpoint the limit where the confinement- induced 
ordering (with zero configurational entropy) gives way to disordered packings (with 
nonzero configurational entropy) as the number of particles increases. Specifically, we 
want to answer the following questions: How often n 2 disks crystallize in nxn square 
arrays with packing fraction 7r/4, when packed under the above- described protocol? 
How does the transition from square to (mostly) triangular symmetry occur? How do 
the confinement boundary effects fade in the large n limit? 

§ Note that a circle can be considered as a regular polygon with an infinite number of sides. In practice, 
however, we have found that polygons with as few as 100 sides produced results indistinguishable from 
those using (mathematical) circles. 
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Figure 1. (Color online) Probability density function for the packing fraction of n 2 
particles in square containers with hard (red or light grey) and periodic (green or 
dark grey) boundaries, (a)-(g) n 2 =9-81, respectively. The low-coverage peaks for 
the probability density in (a) and (b) have been magnified 6x and 3x, respectively, to 
facilitate their inspection. The additional structure (peaks) in the probability density 
functions are associated to special particle configurations for small (n 2 <49) confined 
systems. For comparison, in (h) we show the PDFs for large systems with n 2 =1024. 



3. Results 



Figure [T] summarizes our findings for the probability density function (PDF) of the 
packing fraction for a series of n 2 number of particles in square containers. In all cases, 
we have built the PDF upon 4000 independent packed configurations (full compression 
runs), after assessing in selected cases that a higher sampling did not change the trends 
seen here. The most prominent feature displayed by PDFs is the number of peaks, with 
up to four distinct high symmetry configurations visible with finite probability. This is 
most apparent for systems with few particles (n 2 <49) in hard- wall containers (red or 
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light grey). It is also important to note the lack of any non-symmetric disordered 
packings with packing fractions between those visible peaks visible. The peak at 
packing fraction $=7r/4~0.78 corresponds to nx» square-array configurations!]]] For 
containers with hard- walls (or hard-boundary conditions, HBCs), such square arrays 
dominate the configurational landscape for up to 25 particles, where the second peak 
observed around $~0.74 also becomes significant. Square-array configurations become 
increasingly less probable as the system size increases; for n 2 =36 only 3.7% of the 
samples pack in such configuration. For n 2 >49, nxn square- array configurations have 
vanishing frequency probabilities, and disorded packings are increasingly dominant. 
After the relaxation of confinement effects which strongly encourage crystallization, 
the mean packing fraction increases with the number of particles as previously observed 
for polydisperse systems [T5j[lB]. The shape of the distribution also changes, increasing 
in width as more random configurations begin to dominate. The occurrence of high 
symmetry configurations other than nxn is notable, even for systems as small as 16 
particles, where 22% of the configurations are different from the most probable. This 
number increases to 62% for 25 particles in a square container with hard walls. 

The HBC features discussed above are the result of confinement, that is the 
combination of finite size plus surface effects. By contrast, due to the lack of surface 
effects, the PBC probability density functions show characteristics that are solely due to 
the system's finite size. A comparison between these two set of results allows therefore 
a separation of the two often conflated effects [15]. 

For example, square-array configurations (with packing fraction vr/4) are almost 
never observed for particles with periodic boundary conditions (PBCs). Instead, 
configurations with patches of square and triangular symmetry are favored for small 
systems with PBCs, leading to a bimodal structure in the probability distributions at 
low n. Similarly to O'Hern and co-workers [HMT9], w ho used a completely different 
model and protocol, we have also observed that even-odd effects are present in the 
shape of the distributions for systems with PBCs. 

For larger systems with HBCs, the PDFs become single peaked. We can extrapolate 
the mean value of the packing fraction distribution $ as a function of the number of 
particles N in the container using [T5"|[T7] 

$(iV) = $buik + aiV- 1 , (1) 

with $buik=0.831 and a=— 2.82. Notice that this value of $buik is in the range of the 
corresponding random-close packing limit of 2D frictionless disks [20] . 

However, the presence of square hard-wall boundaries affects the system 
configuration even when the number of particles is very large, despite the fact that 
the interparticle and wall-particle interactions are only contactlike. Due to the 

|| Our numerical calculations of the packing fraction have an accuracy of ±6% for the containers with 
9 particles. This is characteristic of small systems, with the efficiency of an inflationary algorithm 
increasing steadily with the number of particles. The conclusions of our study are independent and 
thus unaffected by the numerical precision of our calculations 
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Figure 2. (Color online) Selected configurations for 9 [(a)— (d)], 16 [(e)— (h)], 25 [(i)- 
(1)], and 36 [(m)-(p)] particles in a square container with hard walls. For the highly 
symmetric configurations, their sampling frequency is shown. The first column from 
the left contains the configuration with the lowest packing fraction, whereas those of 
the last column correspond to the highest packing ones. Configurations (n) and (o) 
represent two typical 36-particle configurations around the peak of the PDF. Under the 
current protocol, our model did not retrieve high symmetry configurations for n 2 >49. 



incompatibility between the symmetry of long-range crystalline order and the symmetry 
of the container, random arrangements occur with greater probability than observed 
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with PBCs. Confinement, therefore, offers an alternative route to study random packing 
of frictionless monodispersed disks that would otherwise crystallize into a triangular 
lattice. 

It is interestingto observe that under the current protocol, our model retrieved 
the optimal packing^ in all cases, with the exception of n 2 =49. For confined systems 
with few particles (i.e., n 2 <49), the optimal packing coincides with the most probable 
packings (i.e., highest frequency peak of the PDF), whereas for n 2 >49, the optimal 
packing is only found in the tail of the distribution. At the same time, rare-event 
configurations (e.g., the two low-density configurations for 9 particles in Fig. [5^a)-^b)) 
were also found during the sampling of the configuration space. 

Some of the highly symmetric configurations, differing from the most probable, have 
small but significant sampling frequencies, for instance, the configuration in Fig. Efg) 
appears ~6% of the times (corresponding to the second peak from the left in Fig. HJb)). 
An interesting continuation to this work would be on possibilities in targeting particular 
configurations by the combined effect of geometrical confinement with the introduction 
of small yet controlled deviations of the particle shape [22]. This is an important 
question, as experiments in the real world deal, firstly, with particles whose shapes 
are far from the mathematical perfection of a circle (or sphere), and secondly, with 
particles that are by necessity always confined within given boundaries. 

We are thus able to conclude that for monodispersed particles, subject to a 
stochastic driving force (e.g., turbulent flow), the optimal configurations are the most 
probable for small confined systems (n 2 <25), whereas for larger systems, optimally- 
packed structures are found only in the tails of the distributions. These distributions, 
in turn, are increasingly dominated by random-packed configurations as the size of 
the system increases. The existence of high-symmetry configurations (other than 
the nxn square arrays) with non- negligible occurrence probabilities is, therefore, a 
direct consequence of monodispersity and confinement. This is an interesting result 
for applications in devices whose functionality is dependent on the state of order, 
regardless of the low-yield rate for symmetric configurations observed from the present 
protocol [2I6H9]. 

4. Summary and conclusions 

In summary, separating the boundary and finite size effects by sampling a large number 
of packing structures with HBC and PBC, we observe that monodisperse disks crystallize 
only in the nxn square configuration for very small systems, n 2 <36, with the triangular 
crystallization more probable above a confinement limit of n 2 =49 (i.e., random packing 

If By optimal packing we mean the densest arrangement of circles in a square container; a mathematical 
problem of geometric optimization with a venerable history that goes back to the work of Fejes Toth. 
Optimal packings of circles in square containers, for a number of circles in the range discussed in 
this paper, have been found by strict mathematical arguments or numerical simulations. An excellent 
review can be found in Ref. [21] 



Confinement effects on monodisperse disk packings 



8 



can be seen as a protoform for triangular lattice formation). Boundary effects, while 
diminishing as the number of particles increases, are still visible even at large n due to 
the mismatch in symmetry between the particle and container. As strict confinement 
has been seen to produce packing configurations not found often in nature, confinement 
can easily be used as a parameter to tune the experimental realization of two dimensional 
structures within colloidal suspensions, microfluidic channels, or molecular superlattices. 

We believe that our results for square containers are indicative of a general behavior 
of confined granular disks and that the main qualitative aspects will be observed as well 
in containers with different shapes. A more difficult and interesting question is whether 
optimal configurations of particles with anisotropic shapes, e.g., ellipses, are also the 
most probable packings under a 2D stochastic agitation protocol. This work also leads 
to several other questions, including whether a packing protocol can be tuned to increase 
the yield rate of nontrivial high-symmetry configurations while staying amenable to 
experimental realization. This includes the essential question of the threshold value 
above which the polydispersity mitigates the ordering provided by confinement. These 
questions surely warrants further investigation, and we shall pursue them in future work. 
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